Symmetry of the charge density wave in cuprates 
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We derive and analyze an eflective Ginzburg-Landau (GL) functional for a charge density wave 
(CDW) for a model of electrons on a tight binding square lattice with density-density interactions. 
We show, using realistic electronic dispersions for the cuprates, that for the simplest GL theory, the 
preferred symmetry is of the uni-directional (stripe) type, but upon inclusion of third-order terms, 
the CDW may exhibit a transition from a uni-directional to a checkerboard pattern depending on 
the strength and range of the interaction. 

PACS numbers: 71.45.Lr, 74.20.De, 74.72.-h 



I. INTRODUCTION 

Modulations of the local density of states in 
cuprates, which were observed by scanning tunnel- 
ing microscopy in Bi2Sr2CaCu2 08-i-5 (Bi2212)i"— and 
Ca2_2;Naa;Cu02Cl2^i^, have attracted significant atten- 
tion. Although a number of theories for these ob- 
servations have been proposed, ^"^^ the nature of the 
modulated state is still debated. The modulations are 
strongest in the underdoped region of the phase diagram, 
a faithful description of which could be a difficult task. 
Fluctuations of the superconducting order parameter, in- 
trinsic disorder, and competing/coexisting order, all po- 
tentially play a role in underdoped samples. 

Nevertheless, short of knowing the exact Hamiltonian 
governing the low energy phenomenology of the cuprates, 
it is instructive to isolate the role of separate contribu- 
tions by focusing on one of them at a time. The goal 
of this article is to study the influence of the electronic 
dispersion in determining the nature of potential charge 
density waves in cuprates. This is realized by deriving 
a Ginzburg-Landau (GL) free energy and analyzing the 
symmetry of the possible charge modulations, with the 
coefficients of the free energy determined from disper- 
sions extracted from angle-resolved photoemission data. 



II. GINZBURG-LANDAU FREE ENERGY FOR 
A CDW 

Our starting point is the generalized extended Hub- 
bard model of interacting electrons on a tight binding 
square lattice. In terms of electron creation operators 
'^k Hamiltonian of the model is: 
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FIG. 1. Diagrams corresponding to the GL coefficients Sa, 
Sb, Sc, and Sd- 



U-V (extended Hubbard) model when g(r— r') is also non 
zero for the nearest neighbor sites r' = r ± ax, r ± ay. As 
a function of the charge-density-wave order parameter 
A(Q), the effective GL free energy (per lattice site) to 
quartic order can be written as 

F(A) ~ F(0) = E - ^(q' 
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H = - A^)^I>k. + 2 E r')n(r)n(r') (1) 



where ^ denotes the energy dispersion (with fi the chem- 
ical potential), r and r' are the sites of the lattice, and 
n(r) is the charge density. The model reduces to the Hub- 
bard model when g{r — r') — goS^.r', and to the so-called 



^ E ^91^92^93^94'^ ( E*M '^4(gi,g2,(?3,g4) 

91,92,93,94 \n=l / 



(2) 



where the summations are performed over the momenta 
q„ and the bosonic Matsubara frequencies fin = 27rn//3, 
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F{Q) is the free energy of the normal state (A^ = 0), 
and the notation = (q„,r2„) is used for brevity. The 
coefficients of the expansion can be expressed through 
Greens functions Gk,w„ = {i^^n ^ + ^J■)^^ as 



where Q _L Q, 



S'4(qi,ili) = GkGk+qiGk+qi+q2 Gk+qi+q2+q3 



(3) 
(4) 
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where k = (k, a;„). After summation over fermionic fre- 
quencies uin = ■7r(2n + l)//3, one obtains the usual result 
for the (2D) charge susceptibility 



a'dk /(^k) - /(^k+q) 
(27r)2 ^k - a+q + 



(6) 



with the integral over the first Brillouin zone. In our 
mean field theory, we only need the cubic and quartic 
coefficient functions 5*3 and 5*4 for — 



a'^dk 



/(fk) 



(■fk - ^k+pi)($k - Ck+pi+pa) 



c.p. 



(7) 

where c.p. denotes the two other terms obtained by cyclic 
permutation of the momenta (k, k + pi,k + pi + P2). 
Similarly, for ^4 one finds 



a^dk 
(2^ 



/(Ck)(^k - ■fk+pil 



(^k - Ck+pi+p2) ^(Ck — Ck+pi- 
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where c.p. denotes the three other terms obtained by 
cyclic permutation of the momenta (k, k + pi, k + pi + 
P2, k + pi +P2 +P3)- As the temperature or the effective 
interaction ^(r — r') varies, for a specific set of equivalent 
momenta Q.^, the quadratic coefficient of the GL expan- 
sion may change sign, and a phase transition takes place. 
In the vicinity of such a transition, it is sufficient to con- 
sider only the coefficients of the GL free energy evaluated 
at momenta Qi and their harmonics. We will be specifi- 
cally interested in situations where these equivalent mo- 
menta lie along the symmetry lines of the Brillouin zone, 
and thus the total number of them is restricted to four 
for the square lattice (typically, the maxima of x will 
lie along such lines). Consider first a simplified expan- 
sion where only the fundamental harmonics at Qi are 
retained. In such a case, the GL expansion reduces to 

F{A)-F{0) = rc,{\Aci\^ + \A^f)+ 
1{2Sa + 45b) ( 27|AQnAQp + |Aq|4 + \A^A (9) 



2x(Q), 



(10) 



fii = is implied, and the coefficients Sa,b,c,d, which 
are described by the diagrams in Fig. [l] are defined as 



5a = 54(Q,-Q,Q,-Q) 
Sb = 54(Q,Q,-Q,-Q) 

5i5 = 54(Q,Q,-Q,-Q) 



(11) 



with 7 = 2Sa+4S° • Fo'^nially identical GL expansions 
have been studiecf in the past, and for a momentum- 
independent interaction have been recently employed by 
Yao et al. in connection with CDWs that occur in the 
rare earth tri-tellurides (Ref. 20). The stability require- 
ment dictates that 25^^ + 45b > and 7 > — 1. Pro- 
vided these are satisfied, the symmetry of the CDW in 
the ordered phase is deterniined^^ by 7: when 7 > 1, the 
free energy is minimized by choosing either Aq = or 
Aq = resulting in one-dimensional stripes. For 7 < 1, 
the minimum of the GL free energy is achieved by choos- 
ing |Aq| — |Aq|, and the CDW has a two-dimensional 
"checkerboard" pattern. 



III. CHARGE SUSCEPTIBILITY AND 
INSTABILITY WAVE VECTORS IN CUPRATES 

The calculations were performed for three sets of six- 
parameter tight-binding fits to the dispersion, based on 
angle-resolved photoemission data, that were previously 
used to model the spin susceptibility^!. The coefficients 
of these dispersions (tbl, tb2, and tb4) can be found in 
that work, tbl is based on earlier photoemission data for 
Bi2213^ and is characterized by a van Hove singularity 
at (tt, 0) that is 34 meV below the Fermi energy, result- 
ing in a sizable anisotropy of the Fermi velocity around 
the Fermi surface. tb2 is based on more recent photoe- 
mission data for Bi2212^ which indicates an isotropic 
Fermi velocity. tb4 is based on photoemission data for 
underdoped La2-2;Sra;Cu042^. 

In Fig. [21 we show x(li ^ — 0)- There are two sets 
of maxima. One set forms a box-like structure around 
the zone center, and it is this set which will be assumed 
to dominate the charge response since we anticipate that 
(7q will be maximal at g = 0. The second set is concen- 
trated around the (tt, tt) point, and is thought to dom- 
inate the spin response since the superexchange inter- 
action is maximal at this point. For the first set, the 
absolute maximum of x typically occurs along the zone 
diagonal {q^ — qy). This can be qualitatively understood 
by employing an argument due to Schulz^. Although 
a bond centered vector does a better job of nesting the 
antinodal sections of the Fermi surface near (tt, 0), a diag- 
onal vector connects twice as many Fermi surface faces. 
Therefore, were it not for the momentum dependence of 
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FIG. 2. (Color online) x{^) calculated from the tight binding dispersions tbl, tb2, and tb4 (left to right}^. Top row: Two- 
dimensional map of x(ci) for T — O.OOlt, where t is the near neighbor hopping parameter of the tight binding fit for ^k- Since 
x(q) is four-fold symmetric, only one quarter of the Brillouin zone is shown. Bottom row: xi'i) calculated for T — 1 meV, 5 
meV, and 10 meV along the various symmetry lines of the zone, with F — (0, 0), X — (1, 0), and Al — (1, 1) in vr/a units. 



the interaction, one would conclude that the CDW typ- 
ically cannot be directed along the bond directions, at 
least within the model considered here. 

In two dimensions, the susceptibility does not diverge, 
even at T = 0, unless the faces of the Fermi surfaces are 
perfectly nested. This is never the case for a realistic 
dispersion, although the Fermi surface does contain ap- 
proximately nested sections that are responsible for the 
sharp peaks in x- This is illustrated in Fig.[3]for the peak 
lying along T — X. The discontinuities of the slope at mo- 
menta qi and q2 (which define the plateau in the insets) 
are due to the points of the Fermi surface k and k' sep- 
arated by such that the tangents to the Fermi surface 
at k and k' are parallel. The lack of divergence in the 
susceptibility implies that the charge-density wave can 
develop only if the interaction strength exceeds a critical 
value. 

The ordering vector in cuprates observed by scanning 
tunneling microscopy lies along the Cu-0 bonds, in con- 
trast with our results for x- On the other hand, the 
true ordering vector is determined by the product of x 
and g. We anticipate that has a maximum at q — 
and falls off smoothly with increasing q. If this fall off 
is relatively isotropic and steep enough, then, as demon- 
strated in Fig. 131 this can lead to ordering along the bond 
direction instead, since q along the diagonal of the box 
structure shown in Fig. [5] is \/2 times larger than q along 
the bond direction. In general, the momentum q = Q at 
which the charge-density instability occurs is determined 
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FIG. 3. (Color online) Susceptibility for T = along T - X 
for dispersions tbl (top left), tb2 (bottom left) and tb4 (bot- 
tom right). In all cases, x is finite due to imperfect nesting. 
The Fermi surface for tbl is shown in the top right panel. 
For dispersions tbl and tb4, the susceptibility has two cusps, 
corresponding to momenta connecting points in a quadrant 
of the zone (shown for tbl by the dots) with their partners 
in the adjacent quadrant. The characteristic bending of the 
Fermi surface near the antinodal direction, which is similar 
for the tbl and tb4 fits, is absent for the tb2 dispersion where 
the susceptibility has a single cusp. In all cases, on one side 
of each cusp at Qxo, the slope of x{Qx,0) is finite, while on 
the other side, x{QxO,0) - xiQ^^O) oc \Qx - Qxo\^^'^ and the 
slope is infinite. 



4 



by the highest temperature such that the condition 

-^=x{ci,T,) (12) 

is satisfied. For concreteness, we assume that the inter- 
action has the following simplified form: 

--^«a + /3|aq|2 (13) 

at small momenta. 

Given this form of the interaction, the first task is to 
identify the momentum Q and the temperature Tc at 
which the CDW first develops. As shown in Fig. [SJ the 
solutions generally fall in four different classes. For a 
fixed curvature /3, there are no solutions provided that 
a is large enough. This was to be expected, as the log- 
arithmic divergence of the susceptibility is cut off due 
to imperfect nesting. Thus, for sufficiently weak interac- 
tions, no ordering occurs even at T = 0. In the opposite 
limit of strong interactions and medium- to long-range 
interactions, the instability is at 9 = 0. This is due to 
the strong reduction in the momentum dependence of x 
as the temperature is raised, as can be appreciated from 
Fig. m In this high T limit, the susceptibility maxima 
are confined to the region near the M point of the zone, 
which leads to solutions for these wave vectors for small 
values of the curvature /3 < 0.1. This is not shown in 
Fig. O The reason is that the parabolic approximation 
we apply for l/gq is only valid near q = Q. We anticipate 
that for larger values of g, the charge interaction is sup- 
pressed, and the magnetic interactions become dominant 
instead. 

We therefore focus on the more pertinent cases of the 
two regions that are shown in Fig. [S] by the green circles 
and blue diamonds. In these cases, the instability occurs 
either at a momentum along V — X (green circles) or 
T — M (blue diamonds); two representative examples are 
shown in Fig. In either case, the ordering momentum 
lies in the immediate vicinity of the "box" structure that 
surrounds the F-point (Fig. [2]). It is straightforward to 
show that the boundary separating the regions shown by 
the green circles and the gray signs in the parameter 
space must be a straight line. This boundary corresponds 
to the limit Tc 0; in Fig. |4] it represents a family of 
parabolas with different a and /3, all passing through 
the maximum of the zero-temperature susceptibility at a 
momentum Qo ~ 0.3557r/a (Fig. [31 tbl). The boundary 
is thus described by a straight line; 

a + pQl = x{T^Q,Qo) (14) 

The other main features of the instability diagram can 
be understood in a similar fashion. 



IV. QUARTIC TERMS AND THE SYMMETRY 
OF THE CDW IN CUPRATES 

Once the ordering momentum Q is known, the coef- 
ficients of the Ginzburg-Landau free energy, evaluated 




M r X 

Wave vector q 

FIG. 4. (Color online) Depending on the momentum depen- 
dence of the interaction gq, a CDW with modulations along 
the Cu-0 bonds (F — X) can be stabilized despite the fact 
that the maximum of the susceptibility along this direction is 
smaller than the peak value along T — M. Thin solid lines de- 
scribe the momentum dependence of susceptibility for the tbl 
dispersion along the contour M — T — X, for temperatures from 
T = 0.001 eV to T = 0.055 eV, with increments of T = 0.002 
eV. The inverse interaction strength — l/2(;q = a-|-/?(aq)^ for 
two different sets of parameters (a, /?) is shown by solid green 
(a = 0.40 eV"\ /3 = 2.38 eV"^) and dashed blue (a = 0.15 
eV-\ P = 2.75 eV"^) parabolas. 

at the transition temperature T — Tc, can be computed 
rather easily. Fig. [S] illustrates the temperature depen- 
dence of the quartic coefficients calculated for the mo- 
mentum Q = Qox, where Qo ~ 0.35 77r/a is the instabil- 
ity momentum for a representative case of a = 2.4 eV~^ 
and j3 = 0.38 eV"^ At high temperatures {T :$> t), the 
leading-order result for all coefficients Si {i = A, B,C, D) 
has the same functional form 48T~^-|-. . ., as can be shown 
by a straightforward expansion of the integrand in Eq. @ 
in powers of £,k/T. In the opposite limit of low tempera- 
tures, the coefficients Sa, Sb, and Sc diverge as T"'^/^, 
while Sb appears to remain finite. At high temperatures, 
7 — 2, independent of any particular form of the disper- 
sion. As can be seen from Fig. [51 7 > 1 is satisfied at 
all temperatures. Therefore, in the simplified description 
where the additional harmonics are neglected, one would 
conclude that ID stripe ordering at a momentum Qq is 
the preferred state. The effect of the harmonics will be 
discussed next. 

Notice that while the coefficients of the quartic term in 
Eq. ^ are nominally determined only by the dispersion 

(Eqs. ^ and ((TT|) ). they also indirectly depend on the 
interaction, ^q, since this enters into the determination 
of the ordering vector as discussed above. However, the 
functional form of the interaction affects quartic coeffi- 
cients in another fashion, which turns out to be rather 
pronounced. When deriving F{A), we left out of the fi- 
nal expression Eq. ([SJ all momenta but those where the 
instability first develops. It happens that the inclusion 
of the third-order terms in the effective action involves 
harmonics of the form A2Q and A^Q_j_Q. These terms 



FIG. 5. (Color online) Instability diagram for the tbl disper- 
sion. The region marked by gray "+" symbols corresponds 
to interactions so weak that even at T = 0, no CDW order- 
ing occurs. For (a, (5) in the parts of the diagram marked 
by green circles (blue diamonds), the instability happens at a 
momentum lying along T — X (F — M) . Even for moderately 
localized interactions (/? small) the instability first appears 
at q = 0, provided that the interaction is sufficiently strong 
(a < 2 eV-M. 
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FIG. 6. (Color online) Temperature dependence of the quartic 
coefficients Sa, Sb, Sc, and So for Q f» (0.3577r/a, 0) (tbl 
dispersion). The inset shows the linear combinations of the 
quartic coefficients that determine the symmetry of the charge 
density wave. The condition 7 = {8Sc + 4:Sd) /(2Sa+4:Sb) > 
1, which is fulfilled for all temperatures in this case, corre- 
sponds to one-dimensional stripes. 



lead to a renormalization of the fourth order coefheient 
and, for a wide range of parameters, this can affect the 
preferred symmetry of the CDW, as shown in Ref. for 
the case of the rare earth tri-tellurides. That is, inclusion 
of these harmonics leads to an additional contribution to 
F(Aq, Aq) of the form 



FIG. 7. (Color online) Left panel: Temperature depen- 
dence of the susceptibility at momenta Q, 2Q, and Q + Q 
for momentum Q ~ (0.3577r/a, 0) (tbl dispersion). Right 
panel: Temperature dependence of the third-order coefficients 
S3(Q,Q,-2Q) and S'3(Q, Q, -Q - Q). 



r2Q(|A2Q|' + |A 

+ (AqA_2Q + A^A 



^2Q 



lA 



-2Q 



+ c.c. 



CQ ( A_Q_Q Aq Aq + Aq_qAqA_q + c.c. 



where 



^'Q = -253(Q,Q,-2Q) 
CQ = -453(Q,Q,-Q-Q) 



1^) 

(15) 

(16) 
(17) 



In the expression for SF3 above, we omitted all terms 
higher than cubic order. These harmonics can be inte- 
grated out, as was shown by Yao et aZ.^", and the result- 
ing correction reads 



bo 

<5F3 = 2.(|Aq| 

r2Q 



|AqI^)- 



lAc 



'Q+Q 



|AqP (18) 



These terms depend on the specific form of the interac- 
tion (7q, which affects the coefficients rq and the choice 
of momentum Q. The coefficients S3 are shown in Fig.[7l 

One can estimate the effect of the corrections by notic- 
ing that the GL coefficients are modified as follows: 



2Sa + 45b ^ 2Sa + 4Sb 



8Sc + 4Sd 8Sc + 'iSo 



453(Q,Q,-2Q)2 



1 



(19) 



+ X2Q 



252Q 

1653(Q,Q,-Q-Q)2 



1 



2.9, 



Q+Q 



^Q+Q 



(20) 
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In the situation shown in Fig. |4l the additional term in 
Eq. (IT51) can be neglected (i.e., the response at 2 Q is typ- 
ically small). On the other hand, the one in Eq. (f20| can 
be significant, given the box-like structure of x about 
the zone center (i.e., the response at Q -f Q can be 
large). As a result, the net effect of the cubic corrections 
on the GL free energy is a reduction of the coefficient 
7 describing the relative magnitude of the mixed term 
lAqplAqp. Consequently, the range of parameters for 
which the checkerboard symmetry is realized could be 
generally increased. 

Evaluation of the renormalized quartic coefficients 
from Eqs. (|19p and (|20l) leads to a refined version of the 
instability diagram shown in Fig. [51 In the region of the 
diagram that corresponds to low transition temperatures, 
the contribution from the quartic coefficients dominates 
over the third-order terms, and the preferable state is 
one- dimensional . 

For smaller curvature /3, as one approaches the part of 
the diagram where the instability occurs along a diago- 
nal wave vector, the contribution of the third-order terms 
becomes more pronounced and results in a transition to a 
two-dimensional state, as described above. Both the one- 
dimensional and the two-dimensional cases are character- 
ized by the dominant order parameter Aq cx {Tc — T)^/'^. 
In all cases, there are subdominant higher harmonics 
of the order parameter at wave vectors 2Q, Q -I- Q, 
etc. that in the vicinity of Tc behave as integer pow- 
ers of Aq. For example, the most pronounced subdom- 
inant order parameter for the "checkerboard" state is 
^Q+Q « IAqIV'^q+q « (Te-r)/rQ+Q. Note that as /3 
is reduced, ?'q+q also decreases (Fig. [J]). Thus for fixed 
Tc — T, the subdominant order Aq^g will be growing, 
and one expects that the simple GL scenario considered 
so far will break down at some point. 

Indeed, for even lower values of /?, the renormalized 
value of 7 falls into the 7 < — 1 range. In other words, 
the third-order terms destabilize our reduced GL free en- 
ergy expression that included only terms up to the fourth 
order in Aq. This instability implies the necessity to in- 
clude higher-order terms in the GL free energy expression 
for that part of the parameter space. These terms would 
generally restore the stability of the GL free energy, but 
would result in a first-order transition that cannot be de- 
scribed in any quantitative way using our approach. At 
the transition temperature Tc, a finite amplitude of Aq 
develops, and unless a weakly first order transition oc- 
curs, one in principle has to include an infinite number 
of terms in the GL free energy to describe it. The cou- 
pling of the order parameter Aq to the modes at 2Q, 
Q + Q, and higher harmonics in this case implies that 
A2Q, Aq^Q, etc. also acquire finite values just below Tc- 
Since none of the higher harmonics have been observed 
in experiment, a detailed analysis of this phase is left for 
future study. 
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FIG. 8. (Color online) A refined version of the instability 
diagram of Fig. [5] Instability at momenta lying along the 
Cu-O bonds can result in a ID stripe phase (green crosses) 
or a 2D "checkerboard" phase (red squares). In the region 
marked by magenta circles, the cubic corrections that involve 
modes with momenta Q+Q destabilize the GL free energy (at 
the quartic level) by violating the stability condition 7 > — 1, 
corresponding to a first-order transition to a state with finite 
Aq, Aq-, Aq^q, and higher harmonics. 



V. CONCLUSIONS 

One faces two difficulties when attempting to recon- 
cile the atomic-scale modulations observed in real space 
by scanning tunneling microscopy (STM) with the en- 
ergy dispersions obtained from photoemission. First, one 
must account for the fact that the modulation wave vec- 
tor observed by STM is directed along the Cu-O bonds, 
while the charge susceptibility extracted from the pho- 
toemission dispersions is largest along the diagonal of 
the zone. At the level of a weak-coupling theory with an 
effective electron-electron interaction, this apparent con- 
tradiction implies that the interaction must be momen- 
tum dependent. Moreover, one obtains rather stringent 
constraints on the range and strength of this effective 
interaction. 

Even when the interaction is such that the instability 
momentum is oriented correctly, there remains the ques- 
tion of a single Q (stripe) versus a double Q (checker- 
board) state. Here, at the level of Eq. (9), we generally 
predict a single Q state. This is consistent with neutron 
and most STM data. On the other hand, corrections 
to this theory due to third-order terms (Eq. (20)) can 
give rise instead to a double Q state for larger values of 
the interaction strength or interaction range /3~^/^ 
(Fig. 8). Interestingly, the well known difficulty of dif- 
ferentiating whether the observed modulation patterns 
arc uni-directional or checkerboard-like implies that al- 
though uni-directional behavior is the most likely, the 
real materials are near the 1D-2D phase boundary^* in 
Fig. 8. 
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